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ABSTRACT 

We present the Hi data for 5 spiral galaxies that, along with their Ha rotation curves, 
are used to derive the distribution of dark matter within these objects. A new method 
for extracting rotation curves from Hi data cubes is presented; this takes into account 
the existence of a warp and minimises projection effects. The rotation curves obtained 
are tested by taking them as input to construct model data cubes that are compared to 
the observed ones: the agreement is excellent. On the contrary, the model data cubes 
built using rotation curves obtained with standard methods, such as the first-moment 
analysis, fail the test. The Hi rotation curves agree well with the Ha data, where they 
coexist. Moreover, the combined Ha + Hi rotation curves are smooth, symmetric and 
extended to large radii. 

The rotation curves are decomposed into stellar, gaseous and dark matter contri- 
butions and the inferred density distribution is compared to various mass distributions: 
dark haloes with a central density core, A Cold Dark Matter (ACDM) haloes (NFW, 
Moore profiles). Hi scaling and MOND. The observations point to haloes with con- 
stant density cores of size rcore fopt and central densities scaling approximately as 

Po c< Tcore ■ ACDM modcls (which predict a central cusp in the density profile) are in 
clear conflict with the data. Hi scaling and MOND cannot account for the observed 
kinematics: we find some counter-examples. 

Key words: galaxies: methods: data analysis - kinematics and dynamics ~ galaxies: 
spiral - dark matter. 



1 INTRODUCTION 

One of the tantalising and yet unanswered questions of con- 
temporary astrophysics is that concerning the nature of dark 
matter (DM). While its existence has been inferred for sev- 
eral decades, only in the recent past its amount has been 
carefully assessed on various astrophysical scales. On galaxy 
scales, the classical and still most powerful tool to unveil 
DM is that of rotation curves. Compare d to the pioneerin g 
work of lRubin. Thonnard fc FordI (Il980l) and lBosmal Jl98ll) , 
modern studies take advantage from greatly improved obser- 
vational techniques, facilitating much more significant vali- 
dations of cosmological theories, such as Cold Dark Matter 
(CDM). Although the prime astrophysical importance of the 
juxtaposition of observati ons and models was as sessing the 
amount of DM in galaxies JPersic fc Saluccil988l) . the timed 
question has become its distribution. This issue directly per- 
tains to the nature of this component, since a precise knowl- 
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edge of the density profile p{R) is likely to be decisive: the 
DM nature is likely to determine the halo properties. 

Yet, a wealth of observations have not been able 
to settle but the most fundamental issues, such as the 
validity of the CDM scenario on gala ctic scales (e.g. 
iNavarro et all Il996l. iMoore et all ll998D . Some of the 
shortcomings on the observational side were related to 
the limitations of optical rotation curves (internal rms, 
limited radial extent, dust contamination), while radio 
studies were affected by the well-known limited spatial 
resolution. Nevertheless, progress has been made, due 
to the increasing number of investigations on dwarf spi- 
rals and low-surface brightness (LSB) galaxies. These 
studies cast serious doubts on one of the fundamental 
properties of CDM haloes, namely their 'cuspiness' in the 
galaxy centres, and infer the presence of large core radii 
in the DM densit y profiles (see e.g. [ Flores fc Primack 
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However, due to the many steps in the data analysis, there 
can be subtle systematic errors that could distort the results 
or in any case render the results very poorly constrained. 
This has triggered off the recent debate concerning the 
reliability of the data and how well the mass models are 
really constrained. There are claims that the observations 
could actually be consistent with the dark matter density 
profiles predicted by the CDM simulations, not only by 
con sidering the Hi data alone l|van den Bosch et al.l bood 
and Ivan den B osch & Swatcrs'200i), but also by combinin g 
Hq and Hi data iPrimack ,2002 and iS waters et alj|2003all . 
This is the reason why particular care should be taken 
in choosing a suited sample and in perfor ming the data 
analy sis. Note that recent simulations (e.g. iNavarro et alJ 
|2003^ do not converge to a well-defined value of the inner 
slope down to the resolution limit (about 1 kpc), even 
though the slope of the dark matter density profile (defined 
as — dlnp/dlnr) at 1% of the virial radius is still about 
1.2 for a typical galaxy. Notice also that the observational 
results on spiral galaxies show a discrepancy with the 
standard ACDM predictions well beyond 1 kpc, i.e. well 
beyond the resolution limit of the simulations. 

In this paper we study a sample of galaxy rotation 
curves ideal for deriving the properties of the dark matter 
haloes around galaxies; it consists of five late- type bulge- less 
normal (high s urface b rightness) spiral galaxies; compared 
to the sample of lBorriell o & Salucci ( 2001) we have the same 
resolution in the central regions, but a larger spatial exten- 
sion allowing us to determine the size of the core radii. One 
of the main concerns is that the rotation curve has both a 
high spatial resolution and a large extension, i.e. beyond the 
optical radius; this is typically achie ved by combining optical 
(H q) and radio (Hi) dat a (see e.g. Ide Blok fc Bosmall2002l 
and lSwaters et a'lll2003al) . The former provide the necessary 
high resolution (1" . . . 2") while the latter allow us to trace 
the potential out to large radii, typically 2-3 times the op- 
tical radius. It has been known for a long time that Hi data 
usually suffer from the lack of resolution, which can affect 
the results, especiall y in the innermost parts; th is problem 
is nicely reviewed bv lvan den Bosch et alJ ll2000l) . 

The difficulty in deriving reliable rotation curves from 
Hi data cubes resides in the fact that a data cube is 4- 
dimensional (RA, dec, flux density and radial velocity), 
while the rotation curve is 2-dimensional (rotation veloc- 
ity vs. radius) . This means that starting from the data cube 
we have to take into account two types of considerations: 
one concerning the positions in which the rotation curve 
should be traced, i.e. how to reduce the spatial dimensions 
from (RA, Dec) to galactocentric radius. This is usually 
achieved by taking the position-velocity diagram along the 
major axis or by fltting concentric ellipses around the cen- 
tre (tilted-ring modelling of the velocity field). The other 
consideration concerns the rotation velocities that should be 
attributed to these galactocentric radii: at the positions de- 
fined above we want to associate only one of the possible 
velocities that could in principle be derived from the spectra 
at single points (flux density vs. radial velocity), i.e. the ra- 
dial component of the rotation velocity at that radius. This 
is usually done via different methods, most of them assum- 
ing symmetry of the proflles, like the flrst-moment analysis 
(the intensity weighted mean) , or the single Gaussian fltting; 
therefore, when the profiles are not symmetric, i.e. in galax- 



ies with a high inclination and/or a poor to intermediate 
resolution, these methods cannot be applied. 

We have developed a new method for both of these con- 
siderations: a) concerning the positions, we traced the rota- 
tio n curve along the ridge of the warp (this step is discussed 
in tVergani et alj |200l . when the velocity field was not suf- 
ficiently sampled (i.e., for the galaxies ESO 116-G12, ESO 
287-G13, ESO 79-G14 and NGC 7339) and the tilted-ring 
modelling of the velocity field was therefore not possible; b) 
for the velocities we only considered the velocity side of the 
profiles opposite to the systemic velocity, and then corrected 
for the effects that artificially broaden the profile, i.e.: the 
turbulence of the ISM, the instrumental velocity resolution 
and the beam-broadening. The present paper concentrates 
on step b). 

When both steps a) and b) were applied we called 
our method WAMET (WArped Modified Envelope Trac- 
ing method); when only step b) was applied and a ve- 
locity fleld was constructed (i.e., when the velocity fleld 
was sufficiently sampled, the case of NGC 1090) we called 
our method MET (Modifled Envelope Tracing method). 
These approaches were tested by constructing artificial data 
cubes based on geometrical models of the Hi disc that 
were iteratively compared to the observed data cubes. The 
MET/ WAMET method proves to provide a better initial es- 
timate of the rotation curve then more traditional methods, 
such as the first-moment analysis and the single-Gaussian 
fitting. Modelling the data cubes results in a powerful way 
to test and improve the de rived rotation curves, as shown 
e.g. bv lGentile et al.1 ll2003h . 

This work is structured as follows: in Section 2 we 
present the sample and optical data that were used in this 
study, in Section 3 we describe the Hi observations and the 
data reduction and Section 4 shows the data analysis; in 
Section 5 we present the results of data analysis, the new 
method for extracting the rotation curves is shown in Sec- 
tion 6 and the comparison between the observed and the 
modelled data cubes is presented in Section 7. In Section 8 
we introduce the different mass models that were considered; 
the results of the mass decompositions are shown in Section 
9 and in Section 10 we draw the conclusions concerning the 
derivation of the rotation curves and the dark matter haloes 
around the galaxies of our sample. 

Throughout this paper we adopted a value for the Hub- 
ble parameter Ho=75 km s"'^ Mpc~^. 



2 THE SAMPLE 

2.1 The relevance of the present sample 

The present sample is ideal for the study of dark matter 
within galaxies, not only for the characteristics of the galax- 
ies, but also for the quality of the data. Table Q presents the 
galaxy sample and its physical characteristics. This sam- 
ple allows us to investigate a range of luminosities approx- 
imately 1.5 mag fainter than L*, ideal for investigating the 
dark matter component, because it starts dominating the 
kinematics at relatively small radii. 

In Table 121 the crucial characteristics of the kinematical 
data are shown: optical data allow us to sample very well the 
region inside the characteristic radii of various dark matter 
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Table 1. Physical parameters of the galaxies. M/ is the I-band absolute magnitude and Mjji is the total Hi mass, rni is the Hi radius 
(defined as the radius at which the surface density drops below 1 Mq pc~^) and is the exponential scale length, r^- is the farthermost 
radius with data and Topt is the optical radius (defined as 3.2rd) and RR us the radio resolution scale. The dynamical mass M-dyn is 
determined at rj. 



Galaxy name 


ESO 116-G12 


ESO 287-G13 


ESO 79-G14 


NGC 1090 


NGC 7339 


Mj (mag) 


-20.0 


-21.7 


-21.4 


-21.8 


-20.6 


Mhi (Mq) 


1.5x10^ 


1.1x10^" 


3.5 xlO^ 


8.5x10^ 


5.5x10^ 


Inclination (°) 


74 


78 


84 


64 


79 




6.7 


7.4 


4.8 


8.8 


4.1 


Tf (kpc) 


11.4 


24.9 


17.7 


29.7 


5.4 


fopt (kpc) 


5.4 


10.5 


12.4 


10.9 


4.9 


r//RR 


5.4 


5.8 


5.0 


10.6 


4.6 


Mdyn (Mq) 


3.3x10^0 


1.9 xlO" 


1.3x10" 


1.8x10" 


3.1x10^0 



Table 2. Optical and radio resolution scales (indicated respectively as OR and RR), and rj, in units of the scale radii of the Burkert 
and NEW density profiles. 



Galaxy 


OR/rcore 


R,R,/'r core 


OR/r, 


RR/ts 


f 1 '^core 


rf/rs 


ES0116-G12 


0.03 


0.49 


0.02 


0.29 


2.7 


1.6 


ESO 287-G13 


0.01 


0.18 


0.02 


0.29 


1.0 


1.7 


ES0 79-G14 


0.04 


0.43 


0.02 


0.24 


2.1 


1.2 


NGC 1090 


0.04 


0.32 


0.03 


0.23 


3.4 


2.5 


NGC 7339 


0.02 


0.15 


0.01 


0.08 


0.6 


0.3 



distributions, while Hi data allow us to probe a region larger 
than these radii, which permits us to investigate the core 
radius phenomenon at large distances and consequently to 
constrain the size of the core radii (note, however, the case 
of NGC 7339 discussed in Appendix A). 

The rotation curves used here consist of the Ha data 
for the inner parts and the Hi data for the outer parts. The 
claim of unreliability of observed k inematics put forward by 
iRhee. Klvpin fc Valenzuelal i2003l) does not apply here: a) 
the derivation of the Hi rotation curve was performed in a 
more thorough and reliable way than the standard tilted- 
ring analysis on the velocity field; b) the galaxies studied 
here were at an inclination optimal to minimise projection 
effects; c) in the regions were they coexist, the Ha and Hi 
data (emerging from different physical processes) do agree, 
implying a high quality of the data^ . The small internal scat- 
ter inside the each radial bin indicates that, for the galaxies 
of our sample, the effects of non-circular motions are neg- 
ligible (even though a more thorough investigation should 
involve 2 -dimensional Ha data ), differently from the case 
shown by'Swaters et al.' ('2003b'). Moreover, it can be shown 
(see e.g. ide Blok. McGaugh & Bosma 2003) that the bias 
towards lower velocity gradients in the rotation curve due 
to a misalignment between the slit with the kinematical axis 
is likely to be small. 



^ Note, however, that this is not entirely true for NGC 1090, in 
which for some radii the Ho lie ~10 km s~^ above the Hi data. 



2.2 Optical data 

The five galaxies of this sample were chosen from the sam- 
ple of 967 galaxies with optical rotation curves presented by 
IPersic fc Saluccil il995ft . according to the following criteria: 
a) they belong to the "excellent" subsample (i.e. the two 
sides are symmetric, the data are extended out to at least 
ropt and the number of data points is ^ 30); b) they have a 
reasonably high total Hi flux; c) they have a relatively low 
I-band luminosity (Mi > —21.8), d) they have a large an- 
gular size, and e) they have an inclination i suitable for HI 
studies (50° < i < 85°). For ea ch galaxy, the raw Ha data 
were binned in groups of 4 to 6 (|P ersic fc Saluccll9'95l) . and 
the given error is the uncertainty on the average value in- 
side the radial bin. We took a "minimum error" equal to 
half the average error in order to avoid data points with un- 
realistically low errors that could bias the rotation curves 
decompositions; for the same reason, when a data point was 
clearly inconsistent with the two neighbouring points and 
the general trend of the rotation curve, the value of its error 
was increased. The complete description of the optical ob- 
servations (spectroscopic and photom etric) can be found in 
iMathewson. Ford fc BuchhornI il992h and iPersic fc Saluccil 
^1995^ . 



3 Hf OBSERVATIONS AND REDUCTION 

We observed ESO 116-G12, ESO 287-G13 and ESO 79-G14 
with the Australia Telescope Compact Array in the 750 m 
and 1.5 km configurations; the resulting baselines range from 
31 m to 1500 m. The galaxies were observed for 12 hrs in each 
configuration. The correlator setup yielded 512 channels of 
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Table 3. Hi Observational parameters 



Galaxy name 


ESO 116-G12 


ESO 287-G13 


ESO 79-G14 


NGC 1090 


NGC 7339 


RA 


3'^13™03^9 


2l'^23'^13=.7 


l''02™51=.2 


2'^46™33=.9 


22''37™46=.9 


Dec 


-57°23'27" 


-45°46'20" 


-65°36'36" 


-0°14'50" 


23°47'12" 


Number of channels 


256 


256 


256 


2 X 64 


2 X 64 


Channel separation (kHz) 


31.2 


31.2 


31.2 


24.4 


24.4 


Synth, beam (arcsec) 


30.4 X 24.4 


32.6 X 24.6 


27.4 X 23.0 


17.6 X 14.7 


13.9 X 13.5 


Synth, beam (kpc) 


2.3 X 1.8 


5.6 X 4.2 


4.0 X 3.4 


3.1 X 2.6 


1.2 X 1.2 


Noise (mjy beam~^) 


1.2 


1.2 


1.1 


0.8 


0.7 



width 15.6 kHz, resulting in a total bandwidth of 8 MHz. 
The data were subsequently Hanning-smoothed to have a 
better signal-to-noise ratio. The final spectral resolution was 
31.2 kHz, corresponding to a velocity resolution of 6.6 km 
s-i. 

NGC 1090 and NGC 7339 were observed with the Very 
Large Array in the C configuration, with baselines ranging 
from 35 m to 3400 m. The galaxies were observed for 8 
hrs with a correlator mode which had 4 IFs, each IF pair 
having slightly different frequency centres; this was done in 
order to find a compromise between velocity resolution and 
bandwidth. We had 2 x 64 channels of width 24.4 kHz and 
a total bandwidth of 2.685 MHz; taking into account the 
overlapping channels, the total number of channels is 112. 
The resulting velocity resolution is ~ 5.2 km s"'^. 

The standard flagging and calibration of the u — v data 
from the ATCA observati ons was performed within the soft- 
ware package MIRIAD JSault et al.l Il99ll) . while the the 
VLA data observations were flagged and calibrated with the 
Astronomical Image Processing Software (AIPS). The rest 
of the reduction was performed in MIRIAD. 

Using a number of line- free channels (10 ... 40) on each 
side of the line, the continuum contribution to the u — v data 
was calculated performing a linear fit; this contribution was 
then subtracted from the data. Then, we pro duced a data 
cube with ROBUST parameters ^lBriggs^ll995^ ranging from 
to 0.3, to reach an ideal compromise between high signal- 
to-noise ratio and small beam size; these values correspond 
to a weighting scheme of the u — v data intermediate be- 
tween uniform and natural weighting. The data were then 
CLEANed, first on a small region with emission and then on 
the whole map, with a total of 2500 iterations per channel 
map. The data were then restored with a Gaussian beam. 
Depending on the galaxies, we obtained beam sizes ranging 
from 13" to 33". The observational parameters are listed in 
Table 01 



4 DATA ANALYSIS 

From the CLEANed cubes, new cubes with lower resolution 
were constructed (with a beam size about twice the origi- 
nal). This was done in order to distinguish genuine emission 
from spurious one; we considered as having reliable emis- 
sion all the regions with an intensity > 2 x cr in the low- 
resolution cube; spurious noise was also removed. Moment- 
maps (total intensity maps) were derived based on the 
"masked" cubes. They are shown in Fig. superimposed 
with optical DSS images. The first contour is the "pseudo 



3-a" l evel defined in a way similar to IVerheiien &: Sancisi 
i200lfl : we define the noise at a certain position in the total 
intensity map as: 

h (N ly 4 H 

where A'' is the number of channels that have been 
added at each position in the total intensity map and 
is the noise in the Hanning-smoothed channel maps. The 
"pseudo 3-ct" level is obtained by averaging all the pixels 
with a signal-to-noise ratio between 2.75 and 3.25. 



5 OBSERVATIONAL RESULTS 

From Fig. we note that the galaxies have a relatively high 
inclination and the Hi emission has a large (up to ~ 3 times) 
extent compared to the optical disc. The Hi masses (Mhi) 
given in Table were calculated from the total intensity 
maps after primary-beam correction using the following re- 
lation: 

Mhi = 2.36 x 10"_D^ J SdV (2) 

where Mhi is in Mq, D is the distance of the galaxy 
in Mpc and J SdV is the total Hi flux density in the to- 
tal intensity map, measured in Jy km s"'^. The Hi masses 
are in agreement with previous single-dish measurements to 
within 20%, implying that we are not losing much extended 
emission due to missing zero spacings. 

The radial distribution of the neutral hydrogen surface 
density is computed by integrating the total intensity map 
over concentric ellipses for NG C 1090 and by means of the 
Lucy method dWarmeljligsl) for the other galaxies. The 
results are shown in Fig. |5| where one can notice that the 
data reach the region where Mhi(< r)/r starts to decrease. 

Since the velocity field of NGC 1090 is sufficiently 
sampled, it cou ld be subjected to tilted-ring modelling 
jBeeemaJliigl) . In this case the MET method was applied 
(see the next Section) to construct the velocity field and 
then the first estimate of the rotation curve was derived by 
performing a tilted-ring fit to the velocity field. In the other 
cases the combination of high inclination and a large beam 
(compared to the size of the galaxy) requires the use of the 
WAMET method: we traced the rotation curve along the 
warp on points defined like in Garcia- Ruiz et ^i^Q 02|) and 
implem ented for studies of the kinematics bv fverganiet alJ 
||200J) on the total intensity map, we fit Gaussians to the 
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Figure 1. Optical DSS images (greyscale) superim posed with the Hi total intensity map (contours). The first contour is the "pseudo 
3-0-" defined similarly to lVerheiien fc Sancis] 120011) and is equal to: 1.0 X lO^^atom cm'^ for ESO 116-G12; 1.1 X lO^^atom cm'^ for 
ESO 287-G13; 1.7 x 10^°atom cm-^ for ESO 79-G14; 1.5 x lO^^atom cm'^ for NGC 1090; 3.4 x lO^'^atom cm'^ for NGC 7339; the 
remaining contours are (15, 30, 45,...) X a, except for NGC 1090 and NGC 7339 where they are (6, 12, 18,...) X a. The beam is shown 
in the upper left corner. 
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Figure 3. A typical spectrum of the galaxy ESO 79-G14 at an intermediate galactocentric distance. The arrows indicate the systemic 
velocity as well as the positions of the velocities derived from the MET/WAMET method, the first-moment analysis and by fitting a 
single Gaussian to the profiles. 



density profiles parallel the minor axis, defining the position 
of tlie peak of tlie Gaussian as the "ridge" of the warp. 

Concerning the analysis of the velocity profiles, partic- 
ular care was taken in deriving the rotation velocity because 
they were in general asymmetric. In Fig. I^we show an ex- 
ample of a typical profile at an intermediate galactocentric 
distance: it is not symmetric and it has a tail towards the 
systemic velocity. This is what we expect in case of a highly 
inclined galaxy, because (see Fig. 01 what we observe is the 
integration along a large portion of the disc, so that material 
with lower radial velocities also contributes to the velocity 
profile. Note, however, that Fig. |^ represents the extreme 
edge-on case: the galaxies studied here were chosen to have 
lower inclinations, so the portion of the disc intercepted by 
the line-of-sight is smaller. 

Moreover, from rece n t results 

iSwaters. Sancisi fc van der HulstJ Il997l iFraternali et alJ 
l200lD . it seems that at least in some spiral galaxies there is 
evidence for neutral hydrogen at several kpc away from the 
galactic plane that rotates more slowly than the gas in the 
disc; this effect also produces an asymmetry of the velocity 
profiles in galaxies that do not have a high inclination. 
The result is that the standard methods (e.g. first-moment 
analysis and single-Gaussian fitting), which implicitly 
assume symmetry of the profiles, provide velocities that 
are biased towards the systemic velocity (see Fig. and 
cannot be used to determine the true rotation velocity in 
the galaxies of our sample. 



6 THE MODIFIED ENVELOPE-TRACING 
METHOD 

In order to solve the problems discussed in the previous Sec- 
tion, we buil t a method similar to the Envelope- Tr acing 
method (e.g. ISancisi fc AllenllToTl ISofue et"allll99'it : the 
aim is to fit only the side of the profiles which we are inter- 
ested in, i.e. the extreme velocity side (with the assumption 
that the gas is present at these extreme velocities). 

The first step was to fit a half-Gaussian from the peak 
of the profiles to their extreme velocity side (i.e. the side 
opposite to the systemic velocity), considering the velocity 
at half maximum, the terminal velocity (Vt)- The rotation 
velocity Vrot at a certain position is then given by: 



Vrot 



\Vt ~ Vsys\ - o.^^isvisMY + [sVobsY + {5Vt,Y 



(3) 



where i is the inclination, Vays is the systemic velocity, 
and the terms indicated with 5V describe the profile broad- 
enings (in terms of full width to half maximum, FWHM): 

5VisM is the broadening due to the turbulence of the 
interstellar medium, by assuming a constant velocity dis- 
persion of the ISM, GisM ~ 10 km s~^ (a typical value for 
spiral galaxies), which implies: 



5Vi 



ISM = OISM' 



/81n2 = 23.5 km s" 



(4) 

Notice that following iKamphui^ (^O^) (i.e. taking 
<yiSM going from 12 to 7 km s~^) would not significantly 
affect the derived rotation curves. 

SVoba is the instrumental contribution to the profile 
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Figure 2. Radial distribution of the neutral hydrogen surface 
density (Shi): the filled triangles denote the approaching side, 
the empty triangles the receding side and the filled circles the 
average. The solid line refers to the ratio Mhi(< r)/r, where 
^'[hi(< ^) is the Hi mass inside radius r. 
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Figure 4. Model of the radial velocities in the plane of the galaxy 
ESO 79-G14 and the regions that would be intercepted by the 
beam in the edge-on case. 



broadening, which we set equal to the channel resolution: 
6.6 km for the galaxies observed with the ATCA and 
5.2 km s~^ for those observed with the VLA. 

5Vh is an estimate of the broadening of the profiles due 
to the beam: as is evident in Fig.^] a larger beam will sample 
a larger portion of the velocity field and will thus broaden 
th e profiles, even on thei r extreme velocity side. According 
to lSancisi fc AUeiil (^23): lower limit to 5Vh can be given 
by setting it to zero, yielding an upper limit to the rotation 
ve locity. To est imate an upper limit to SVb, we assume as 
in iBrauiil l)l997f) . that the expected beam uncertainty in the 
profile widths is 2 • [V{r ± 6^/2) - V{r)], where the "-|-" 
corresponds to a positive gradient of the velocity field and 
"— " to a negative one; Ot is the beam FWHM. 

With an upper limit and a lower limit to SVb, we decided 
to estimate SVb using their average, i.e: 

SVb = V{r±eb/2)-V{r) (5) 

MET/WAMET has the following advantages: 

• The profiles are asymmetric, thus methods assuming 
symmetry cannot be used to derive the rotation velocity; 
they would underestimate the rotation velocities, especially 
in the inner parts. 

• We are only interested in the extreme velocity side of 
the profiles because it is not (or little) affected by projection 
effects. 

• It also estimates the correction for the beam broadening 
of the profiles, which can be substantial in regions where the 
gradient of the velocity field is high. 

The implemented method (WAMET) for warped galax- 
ies can be applied to galaxies where the sampling of the 
velocity field is poor, and it enables us to trace the ro- 
tation curve along any possible warp instead of keeping a 
fixed position angle, like in methods based on the analysis of 
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the position-velocity diagram; the sh ortcominKS of keepin g 
a fixed position angle are discussed in I Vergani et alJ i2004l) . 



7 ROTATION CURVES AND MODELS 

In order to derive the rotation curves, we applied the tilted- 
ring modelling of the velocity field for NGC 1090. For the 
other galaxies we determined the kinematical centre and the 
systemic velocity by minimising the differences between the 
two sides. The errors are the maximum of the three fol- 
lowing values: the difference between the velocities of the 
approaching and the receding sides, our correction for beam 
broadening of the profiles, and a "minimum error" equal to 
{2/ sin i) km s~^. 

The MET/WAMET rotation curve served as an input 
to construct mo del data cubes with the task GALMOD 
within GIPSY Jvan der Hulst et alJ Il992l) : these models 
have been then compared to the observed data cubes. Such 
models are built with the assumption that the neutral hydro- 
gen moves along circular orbits; a series of geometrical and 
physical parameters, among which the rotation curve, allows 
us to create model observations of the Hi discs. These mod- 
els, when compared to the real data cubes, verify whether 
they are a fair representation of the Hi disc. Some input 
parameters are inferred directly from the observations and 
kept fixed (e.g., the surface density profile), while others 
(e.g., the kinematical centre and the systemic velocity) are 
first estimated from the observations and then adjusted in 
order to optimise the match between the observed and the 
model cubes. 

In two cases (ESQ 287-G13 and NGC 1090, the galaxies 
with the steepest rotation curves) the rotation curve derived 
with MET/WAMET had to be slightly modified in order to 
better reproduce the observed data cubes. We found out 
that it was necessary to increase the value of the rotation 
velocity of the 2-3 innermost points by values of the order 
of their 1 cr error; this is probably due to the fact that, 
even with the MET/WAMET method, when the rotation 
curve rises steeply we are not able to retrieve precisely the 
true rotation velocity; this is however possible by modelling 
the data cube. In Fig. |3 we display some channel maps 
(observed and modelled) to show that with our final rotation 
curves we are able to provide a very good reproduction of the 
observations. All the galaxies of our sample were successfully 
modelled with a small warp in the outer parts or with no 
warp at all. In particular, we found necessary to include in 
our model parameters a 5° clockwise warp in ESQ 116-G12, 
a 2° clockwise warp in ESO 79-G14 and a 2° anti-clockwise 
warp in NGC 1090. Concerning the inclination as a function 
of radius, it was kept constant at the values given in Table 
except for the galaxy ESO 79-G14, where a 5° decrease of 
the inclination in the outer parts was necessary to reproduce 
the shape of the extreme channels. The errors on the derived 
orientation angles are of a few degrees l^entilc ct al. 2003). 

Using the rotation curves derived from other methods 
provided data cubes that were clearly different from those 
observed. In Fig. M for the galaxy ESO 116-G12, we show 
some channel maps and the position-velocity diagram along 
the major axis of the observed data cube, of a model created 
from the corrected MET/WAMET rotation curve and of a 
model created from the first-moment rotation curve. The 



latter clearly fails in reproducing some of the observed fea- 
tures, like the shape of the extreme channels and the slope 
of the position- velocity diagram. Fitting single Gaussians to 
the velocity profiles gives very similar results to the first- 
moment rotation curve. 

In order to assess the validity of a model we used a num- 
ber of tools to visualise the different properties of the data 
cube. In particular, the model has to simultaneously match 
the intensity and shape of the emission in the single chan- 
nel maps, the column density distribution in the moment-0 
map, and the intensity and shape of the emission in the 
position-velocity diagrams, along the major and minor axis 
and on slices displaced with respect to the galactic centre. 
The comparison of model data cubes to the observations is a 
very powerful and crucial test to check quality and reliability 
of the rotation curves. 

In Fig. |7| we show the Ha and Hi rotation curves pro- 
jected onto the Hi position- velocity diagrams of the galaxies 
of our sample: the agreement between the two datasets is 
good and, as should be expected, the Hi rotation curves fol- 
low the shape of the last contours of the position-velocity 
diagrams. The Hi rotation curves agree with the Ha data 
within 1 cr for the bulk of the points. In the case of NGC 
1090, however, at some radii there is an offset of about 10 
km s^^ between the two data sets: in Fig. |3 one can see 
that the Ha data lie above the Hi rotation curve; the fig- 
ure shows that it is not an error in the Hi data reduction, 
but rather an intrinsic difference, maybe due to non-circular 
motions associated with the small bar visible in the optical 
image. The optical image of ESO 116-G12 looks a bit asym- 
metric, but the kinematics is not, while the other galaxies 
look symmetric in both the kinematics and the photometry, 
which tells us that apart from NGC 1090 the effects of pos- 
sible hidden bars are negligible. We plot about two points 
per beam for the Hi rotation curve: this was done in order 
to have a comparable number of Hi and Ha points. Plotting 
fewer Hi points would result in the Ha data outnumbering 
the Hi data and biasing the rotation curve fits. 

We also compared the Ha data with the Hi rotation 
curves derived with two different methods (Fig. (HJ: our cor- 
rected MET/WAMET method and the first- moment. The 
latter yields a much worse agreement with the Ha rotation 
curves, especially in the inner parts where the profiles are 
most asymmetric: in these cases the rotation velocities are 
severely underestimated. This shows the importance of the 
particular care that should be taken in deriving rotation 
curves. 

The symmetry and smoothness of the rotation curves in 
the selected sample, along with the good agreement between 
the two datasets make these galaxies ideal laboratories for 
a dark matter study based on the kinematics: this holds 
in particular because they lack large non-circular motions, 
a major factor of disturbance and uncertainty in this kind 
of study. The agreement between the Ha and Hi rotation 
curves also tells us that any systematic errors that might ar- 
tificially render the rotation curve shallower (i.e. extinction 
in the Ha data, bad positioning of the slit, low resolution of 
the Hi data) are small. 

Note that the Ha data are absolutely necessary for a 
mass decomposition such as the one presented in the next 
sections. The uncertainties in the inner parts of the Hi rota- 
tion curves are quite large (typically 10-20 km s"'^), there- 
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Figure 6. Comparison between the observed data cube and two model data cubes for ESQ 116-G12: the MET/WAMET rotation curve 
(left) and the first-moment analysis rotation curve (right); the black contours trace the observed data cube and the grey contours trace 
the model data cubes. On the top we show some representative channel maps (the heliocentric radial velocities are shown in the upper 

left corners) and at the bottom the position-velocity diagrams along the major axis. The contours are -2.5, 2.5, 10, 20, 50 mjy beam"-*^. 
The cross indicates the centre of the galaxy and the arrows show some key positions where the difference between the two model data 
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Figure 8. The rotation curves: the Ha (filled circles) rotation curves, the Hi corrected MET/WAMET rotation curves (filled triangles) 
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curves. 
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fore nearly all the mass models fitted would be consistent 
with the data, disabling us from putting any constraints on 
the dark matter distribution. The actual uncertainties on the 
rotation velocities are probably smaller, due to the precision 
than can be reached with the data cube modelling; however, 
they are difficult to estimate and rather subjective, therefore 
we keep the above estimated errors. 



8 MASS MODELS 

The sample allows us to test in detail different mass mod- 
els with various distributions and relative amounts of dark 
matter. 



8.1 Luminous matter 

In order to study the properties of luminous and dark matter 
in these galaxies, we model the circular velocity V{r) in 
terms of the disc, gas and halo components: 



(6) 



Vdi3c{r) is the contribution of the stellar disc, which 
is scaled according to the chosen (stellar, I-band) M/L ra- 
tio (hereafter T'). In this way, the molecular gas is auto- 
matically taken into account, since it is reasonable to as- 
sume that it is distribut ed like the stellar component (e.g. 
ICorbeUi fc Saluccill2000D . The available photometry shows 
that the luminous matter in these galaxies is distributed like 
in an exponential thin disc, whose va lue is given in Tabled 
l|Mathewson. Ford fc Buchhornll 992\ The exception to this 
is ESO 79-G14, whose profile is significantly different from 
exponential: for this galaxy the actual surface brightness 
profile was used. T\ (assumed to be constant with radius) is 
a free parameter in each model, however, with the condition 
of being larger than 0.2. 

Vgasir) is the contribution of the gaseous disc, derived 
from the Hi surface density distribution, scaled up by a fac- 
tor 1.33 to account for primordial helium; the uncertainty of 
the galaxy distance is irrelevant for the present study since 
at no radius the Hi contribution dominates the kinematics. 

Vhaioir) is the contribution of the dark matter halo, 
derived under a number of assumptions as discussed below. 



8.2 Dark matter 

8.2.1 Burkert halo 

iBurkertj 1^9^ proposed an empirical density distribution 
of dark matter: 



pB(r) = 



PO'^core. 



(7) 



Here, po is the central density and Vcore the core radius. 
This halo has a constant density core, pB^r) — > po for r ^ 0, 
and a decrease of the density as for r — ^ oo. 

The halo has two free parameters: the core radius a = 
Tcore/ropt (i.e., in units of the optical radius, Vopt ~ 3.2rd) 
and the central density po- 

Other cored distri butions were investiga ted (the 
pseudo-isothermal halo, Ivan Alba da et al l Il985l and the 
Universal Rotation Curve halo. lPers'k^ahiccrfcSteilliggd 



hereafter PSS), but they yield results which are almost in- 
distinguishable from the Burkert halo. We include them in 
Table 01 for completeness. 



8.2.2 Navarro, Frenk and White (NFW) halo 

iNavarro et"ai] lll99d) fitted the outcome of N-body CDM 
simulations with a density distribution given by: 



PNPw(r) = 



{r/rs){l + r/rs)^ 



(8) 



Here, ps and are the characteristic density and 
the scale radius of the distribution. They are in princi- 
ple independent, bu t recent results jBuUock e t al. 20(3, 
IWechlser et al.ll2002ll show a correlation between these pa- 
rameters. The NFW density distribution is then a one- 
parameter family, n amely the virial mass M^ir. From 
IWechlser et"ai] ll2002ll we take the relations linking Mvir to 
the concentration parameter c (= rvir/rs), ra and pa, at red- 
shift z = and for a Universe with A = 0.7 and Q,o = 0.3, 
starting with M^uir = |7rA„irPc?'?ir (where A^ir is the virial 
overdensity and its value is about 337 at z = 0, pc is the 
critical density of the Universe and r„ir is the virial radius): 



20 



1O"M0 



: 5.7 



IO11M0 



101 



3 ln(l + c)-^' 



kpc (9) 



(10) 



NFW halo has then a central density cusp, with pnfw oc 

for r ^ 0, and a profile/amplitude which is controlled 
by a free parameter Mvir. 

Notice that in principle, adiabatic contraction of the 
primordial dark matter halo due to baryon infall should be 
taken into account, but since the effect is to render the halo 
even more concentrated, aggravating thus the known prob- 
lems of the NFW haloes, we neglect it. 

We constrain the virial halo mass to be M^ir < 8 x 10^^ 
M0, in that, for a low luminosity spiral, it must presum- 
ably be substantially lower than that of the Milky Way 
and other very luminous gal axies, for which it is safely esti- 
mated : M^ir ~ 2 X 10^^ M0 llChenealur. Salpeter fc TerzianI 
ll99llWilkinson fc Evansll999^ . This constraint affects onlv 
N7339 and ESO 79-G14, due to the relatively Umited exten- 
sion of their Hi rotation curves, which prevents to rule out 
large ACDM haloes. 



8.2.3 Moore halo 



Recent numerical simulations bv lMoor" et al.l (Il998h yielded 
a more concentrated density profile: 



PMooro('") = 



(r/r,)i-5(l-f (r/r,)i'5)' 



(11) 



where pa and ra are the characteristic density and the 
scale radius of the distribution. This density distribution 
has an even steeper cusp (pinoore oc r~^'^ for r — > 0) than 
the previous one. Similarly to the NFW halo, we consider 
this profile as hav ing only one free parameter. Following 
iMoore et al.1 (|l999r) . we define CMooro as being 1.8 times 
smaller than cnfw; it is then derived from Eq. |2l For a 
given virial radius, the scale radius ra of the Moore halo will 
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then be 1.8 times larger than its corresponding quantity for 
the NFW halo, ps can be derived from: 

„3 



101 



2 ln(l + ci-5)' 



(12) 



Also in this case we constrain the virial halo mass to be 
lower than 8 x lO" Mr.^. 



8.2.4 Hi scaling 

Early studies of Hi rotation curves jBosmalllQsih noted the 
fact that the ratio between the Hi surface density and the 
dark matter surface density is approximately constant in the 
outer parts of galaxies (but see Corbcll i fc Salucci 2QO0i) . 
This led to the hypothesis that dark matter could in some 
way be associated with the Hi disc and distributed in the 
same manner; this is what is reasonable to expect in the case 
of models considering f or instance H2 clumps as a compo - 
nent of dark matter SPfenniger. Combes fc Martinetlll994l) . 
In this case the scaling factor for the Hi contribution to the 
rotation curve is a free parameter. 

8.3 MOND 

Accordi ng to MOND, the law of Modified Newtonian Dy- 
namics jMilgroml llQS^. there exists a certain acceleration 
ao below which Newton's law of gravity is no longer valid 
and the expression for the gravitational acceleration reads: 



gir) = [GM{r)ao]'^^/r, 



(13) 



for the 



stellar and 

2 



where M(r) account 
jaseous components and ap = 1.2 x 10~* cm s 
Beeeman. Broeils fc Sander3ll99J) . 



9 MASS MODELLING RESULTS 

The fits were performed by a x'^-minimisation, considering 
both the rotational velocities and their logarithmic gradients 
(V = iSHlUM") which bear a crucial information on the 
matter distribution in a galaxy fsee lPersic fc Saluccilll99CI) . 
The total value to be minimised then is Xtot ~ Xvei +Xv- 
It is worthwhile to point out that the values should only 
be considered as a way to compare the different fits within 
the same galaxy, rather than a probability indicator, because 
the choice of the error bars is quite subjective and we plot 
two points per beam, so the points are not independent; the 
goodness of a particular mass model is also related to the 
fraction of observational points that it hits within 1 cr as 
well as the ones that it badly misses. 

In Figs. IUlto ll3l we show, for each galaxy, the results of 
the fits, the residuals {Vobs — Vmodei) of the fits and the 1 a 
probability contours in parameter space. The case of NGC 
7339 will be discussed in Appendix A. 

The Burkert profile - so as any cored profile - has the 
best fits to the rotation curves, with no systematic devia- 
tion from the observed rotation curves seen in all galaxies. 
None of our ~ 100 data points (considering the five galaxies 
together) is inconsistent with this model, having a residual 
larger than 3 a (where a is the observational error) . The stel- 
lar I-band mass-to-light ratios, which lie between 0.5 and 
1.8, are consistent with population synthesis models (e.g.. 
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Figure 14. Relation between the Burkert central den sity and 
Burkert core radius; the solid line is the relatio n found bylBurkert 
1 1995) and extended to spiral galaxies by ISalucci fc Burkert 
12000): po = 5 X 10~^'^r ~i'l^e-('"'^°'-°/^^)^ g cm-^. The open cir- 
cle represents DDO 47 iSalucci et alJl20n,'ll and the filled circles 
are the galaxies of the present sample. 



iBell fc de JondboOlT) . The core radii are in the range (0.7 - 
2.3) X Vopt, and the central densities are between (0.4 - 3) x 
10"^'* g cm^"'. In Fig. 1141 we plot the gala xies of our sample 
in the po — rcore p lane of |Burkertl ^ll995^. s lightly adapted 
to spiral galaxies bv lSalucci fc Burkerti J200ol) : despite a cer- 
tain scatter, they roughly follow the relation, which certainly 
has an implication for the nature of dark matter. 

The minimum x^ values for the NFW haloes are sig- 
nificantly higher than for the Burkert haloes. The former 
fail to reproduce both the velocities and the shape of the 
observed rotation curves. Moreover, there is a systematic ef- 
fect in the predicted velocities, in the sense that the NFW 
haloes predict velocities in the central parts that are too 
high. A significant number of data points (9) is totally in- 
consistent (i.e., residuals larger than 3 a) with the fits per- 
formed with the NFW halo; if we consider the points with 
residuals larger than 2 a their number increases to 22. An 
additional problem of the NFW haloes is that the emerging 
Tt are generally small and in one case unacceptably low; in 
this case the extreme lower limit of 0.2 is reached. We also 
attempted to reproduce the data by leaving c as a free pa- 
rameter, instead of using Eq.|S| no appreciable improvement 
of the fits is found. Moreover, the values of c and M^ir are 
in some cases completely unrealistic, in particular we obtain 
very low values of c. 

The Moore haloes provide even worse results, with very 
high x^ values compared to the cored haloes, and veloci- 
ties that are far too high in the inner parts, which are a 
consequence of the predicted steep central cusps. The in- 
ability of the Moore haloes to reproduce the observed kine- 
matics appears also in the large number of points (25) hav- 
ing residuals larger than 3 cr. Notice that recent simulations 
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Figure 9. Mass models for the galaxy ESO 116-G12. The solid line represents the best- fitting model, the long-dashed line is the 
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respectively. 1 kpc corresponds to 13'.'4. Below the rotation curves we show the residuals (Vohs ^ ymodel)- At the bottom we plot the 
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Figure 10. Mass models for the galaxy ESO 287-G13 (see Fig. |5]for an explanation of layout and symbols). 1 kpc corresponds to 5'.' 
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Figure 11. Mass models for the galaxy ESO 79-G14 (see Fig.|5]for an explanation of layout and symbols). 1 kpc corresponds to 6'.' 
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Figure 12. Mass models for the galaxy NGC 1090 (see Fig. |5]for an explanation of layout and symbols). 1 kpc corresponds to 5'.'7. 
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Figure 13. Mass models for the galaxy NGC 7339 (see Fig. |5]for an explanation of layout and symbols). 1 kpc corresponds to 11'.' 
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Table 4. Best-fitting results: parameters and associated errors, reduced x^- o, = Tcore/fopt, Po is in units of 10~^* g cm~3 and M„i^ is 
in units of 10^^ Mq. 



Halo typo Parameter ESO 116-G12 ESQ 287-G13 ESO 79-G14 NGC 1090 NGC 7339 



Burkert 
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„ __ + 0.21 
0'™-0.1S 


2-29_o.41 


-,-+0.08 
0-67lo.o7 


0.813q 
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.22 


l-60l~79 
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3.i;^.3 
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„ ..„+0.09 
67 

-o.os 


0-56-0 


31 
.16 
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1 47 
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4 2+2-'' 
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^-''-0.3 





.1 
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1 S5 + 0-1-* 
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1.1 


pseudo-isothermal 


a 
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1.1 
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1.4410-^0 ,.5,+0.34 g20l°-«l 4.8ll°-- 8.00l°-°? 

0.20t°-- 0.56;-- 0.651-- 0.911°;- 0.931°;- 

^ 10.6 6.6 13.0 2.7 5.3 



seal fact 



i.49 



+ 0.28 
0.28 



1.28 
18.8 



+ 0.06 
0.05 



,+0.13 
-0.14 



-0.04 
11.1 



6.76 
2.32 
2.8 



+0.23 
0.23 



+ 0.07 
0.07 



11 63+0-0^ 
ll-t'-^-LOS 



-0.05 
1.1 



t;, 0.811-;- 1.481-;- i.aC,.^^ 



O + 0.02 „.+0.02 „„ + 0.02 „ +0.03 

S-0.02 1-24-0.02 1-23-0.02 2.15 
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jNavarro et al.ll200 3>) predict, at the innermost resolved ra- 
dius, a slope of 1.2 (approximately intermediate between the 
NFW and the Moore profile) in the radial density profile of 
a typical galaxy and that they rule out the Moore profile. 



The Hi-scaling and the MOND models give reason- 
able results in a fe-w cases, but -with some decisive counter- 
examples: we have a significant number of points that are 
inconsistent with the fits in both cases. In particular, the 
present data have a resolution such that MOND is not able 
to fit them. The values of t he Hi scaling factors a re consis- 
tent with previous studies llHoekstra et alJ|2001^ ■ Leaving 
the MOND acceleration parameter ao as a free parameter 
only slightly improves the quality of the fits: they are al- 
ways significantly worse than in case of the cored haloes. 
The same holds if we leave a reasonable (25%) uncertainty 
on the adopted distance. 



10 CONCLUSIONS 

We have presented an analysis of Hi data cubes of a sample 
of five low- luminosity spiral galaxies that are used to derive 
the distribution of dark matter within these galaxies, along 
with available Ha and I-band photometric data. 

We found a new method (the MET/WAMET method) 
for deriving the Hi rotation curves, which was needed to 
take into account a) a position angle that varies with radius 
(this is needed as the tilted-ring modelling of the velocity 
field could be used in only one case), and b) the asymmetry 
of the velocity profiles, because the galaxies of our sample 
have a high inclination, and projection effects yield asym- 
metric velocity profiles. Corrections were applied in order 
to account for any effects that artificially broaden the Hi 
profiles. 

The rotation curves derived in this way passed a crucial 
test, i.e. the construction of model data cubes using these 
rotation curves as an input. In some cases the Hi rotation 
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curves derived from MET/WAMET required a small modi- 
fication, mostly in the inner parts. After these small correc- 
tions the model data cubes are in very good agreement with 
the observations, contrary to model data cubes constructed 
using rotation curves derived from other "classical" meth- 
ods. Our final Hi rotation curves are consistent with the Ha 
data, but up to three times more extended. 

In galaxies well suited for inferring the distribution of 
dark matter, we were able to trace the total potential via 
the rotation velocity and to derive the density distribution 
of dark matter. 

From the analysis of combined Hq and Hi rotation 
curves we reach the following conclusions: 

• Cored haloes uniquely fit the rotation curves, with core 
radii of order ropt and central densities of order 10~^* g 
cm"^. 

• NFW haloes cannot represent the dark matter haloes 
around the galaxies of our sample: they fit badly and a sig- 
nificant number of data points is inconsistent with the fits. 

• The situation is even worse for the Moore model, as 
implied by its steeper inner profile. 

• There are cases in which for radii in the range (1 — 
4)rd the Burkert and NFW halo rotation curves are similar; 
this may explain w hy in the literature cases can be found 
(e.g. Ijimenez. Verd e & Oh 2003} in which both models fit 
equally well a rotation curve of small extension and limited 
resolution. 

We thus reach the conclusion that the galaxies of the 
present sample are uniquely successfully fitted by cored 
haloes, with a core size comparable to the optical radius. 
This suggests the existence of a well-defined scale length in 
dark matter haloes, linked to the luminous matter, which 
is totally unexpected in the framework of CDM theory (see 
ISalucci fc Burkert 200aV 
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APPENDIX A: THE CASE OF NGC 7339 

The galaxy NGC 7339, due to the small extent of the kine- 
matical data and the value of its dark mass, shows a degen- 
eracy between the predicted Burkert and NFW halo rotation 
curves. As can be seen in Fig. ll3l and in Table|lJ both models 
provide equally good fits to the rotation curve. This is not 
due to the fitting procedure but it derives from the fact that 
both models have similar predictions. We treat in detail this 
case to show that ACDM haloes, for certain virial masses 
and in certain ranges of radii has a radial profile that meets 
that of the Universal Rotation Curve (see PSS). 

First, let us notice that from Tabled that the rotation 
curve of this galaxy has a small spatial extent, i.e. about 5 
kpc. Then in Fig. (left panel) we show the NFW halo 
rotation curv e, very different in gene ral from the derived 
actual haloes dSalucci fc Burkerll2000l) . for objects having a 
number of virial masses, ranging from 1 X 10" M0 to 3 X 10^^ 
Mq. However, as shown in the right panel of Fig. lAll for 
virial masses close to about 1 x 10^^ Mq, and for radii up 
to 5 kpc, the NFW and Burkert halo rotation curves almost 
coincide; this similarity, combined with the uncertainty on 
TJ;, is the reason for the equal quality of the two fits. 

With the data used here we cannot probe the regions of 
NGC 7339 where the difference is noticeable. Therefore, for 
galaxies having V{ropt) approximately in the range 120—170 
km s~^ and for the radial range 1 < r/rd < 4, in order 
to distinguish between the two cases, we need kinematical 
information at larger radii. At smaller and larger V{ropt) 
the two profiles are intrinsically much more different and 
the comparison is easier (see PSS). 
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Figure Al. Left: NFW halo rotation curves for galactic haloes; the values of log(M„i^/MQ) are indicated on the plot. Notice the very 
large differences with respect to the Universal Rotation Curve (Fig. 6 of PSS). Right: NFW halo rotation curve (solid line) for an object 
with Myir = 8 X IO^^Mq and a Burkert halo rotation curve with rcore = 6 kpc and po = 4.5 x 10~^* g cm~^, for radii up to 5 kpc. 
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